arXiv: 1503.02095vl [nucl-ex] 6 Mar 2015 


Cosmogenic-neutron activation of Te0 2 and implications for 
neutrinoless double-beta decay experiments 


B. S. Wang, 1 E. B. Norman, 1,2 N. D. Scielzo, 2 A. R. Smith, 3 K. J. Thomas, 1,3 and S. A. Wender 4 

1 Department of Nuclear Engineering, University of California, Berkeley, California 94720, USA 
2 Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
3 Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 
4 Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
(Dated: March 10, 2015) 

Flux-averaged cross sections for cosmogenic-neutron activation of natural tellurium were mea¬ 
sured using a neutron beam containing neutrons of kinetic energies up to ~800 MeV, and having an 
energy spectrum similar to that of cosmic-ray neutrons at sea-level. Analysis of the radioisotopes 
produced reveals that 110m Ag will be a dominant contributor to the cosmogenic-activation back¬ 
ground in experiments searching for neutrinoless double-beta decay of 130 Te, such as CUORE and 
SNO+. An estimate of the cosmogenic-activation background in the CUORE experiment has been 
obtained using the results of this measurement and cross-section measurements of proton activation 
of tellurium. Additionally, the measured cross sections in this work are also compared with results 
from semi-empirical cross-section calculations. 


I. INTRODUCTION 

Neutrinoless double-beta (Ov/3/3) decay [1-3] is a long 
sought-after second-order weak process in which a nu¬ 
cleus (A,Z) transitions to a nucleus (A,Z+2) through the 
emission of two electrons. This process is hypothesized 
to occur only if neutrinos are Majorana particles. Ob¬ 
servation of 0 u/3/3 decay would not only establish that 
neutrinos are Majorana fermions, but may also constrain 
the neutrino-mass scale and hierarchy, and demonstrate 
that total lepton number is not conserved. 

In experiments searching for Oz//3/3 decay, the signa¬ 
ture of interest is a peak at the double-beta decay Q 
value (Q pp). As Oz//3/3 decay would be a rare process, 
minimizing the background rate around Qpp is essential 
for improving the experimental sensitivity. Therefore, a 
detailed characterization of all potential sources of back¬ 
ground is important, as any event that can mimic or ob¬ 
scure the Oz//3/3-decay peak is problematic and must be 
well-understood and, if possible, eliminated. 

To miminize external backgrounds, 0z//3/3-decay exper¬ 
iments operate in underground laboratories, where large 
overburdens decrease the flux of cosmic rays by orders of 
magnitude relative to the flux above ground [4[. Further 
reduction of the remaining cosmic-ray background can 
be achieved with muon-veto detectors, and backgrounds 
from natural radioactivity in the laboratory environment 
can be alleviated with proper shielding. 

Radioactivity present within the detector itself can 
provide a source of background that is difficult to elimi¬ 
nate. Oz^/3/3-decay experiments devote a great deal of ef¬ 
fort into making ultraclean and ultrapure detector mate¬ 
rials free of primordial radioisotopes. However, no matter 
how clean or purely produced the materials are, cosmo- 
genic activation will generate some radioactivity while 
the materials are at or above the Earth’s surface dur¬ 
ing storage, production, or transportation [3HZ|. The 
background contribution from this radioactivity can be 


minimized by ensuring detector materials spend as little 
time above ground as possible and by avoiding air trans¬ 
portation, as the cosmic-ray flux increases significantly at 
higher altitudes EM- At sea-level, activation is primar¬ 
ily caused by the hadronic component of the cosmic-ray 
flux, which is dominated by neutrons m- 

This work investigates the backgrounds associated 
with cosmogenic activation of tellurium, which are impor¬ 
tant to understand for experiments such as the Cryogenic 
Underground Observatory for Rare Events (CUORE) 
ED and the Sudbury Neutrino Observatory Plus (SNO+) 
[T2] that are searching for the C )u/3/3 decay of 130 Te, but 
to date are poorly characterized due to a lack of data. 
As 0z//3/3-decay experiments run for several years, typi¬ 
cally only long-lived cosmogenic isotopes (i.e., that have 
half-lives of order a year or longer) with Q values near or 
greater than the 130 Te Qpp of 2528 keV fT3HT6] will be 
potential sources of background at the 0z//3/3-decay peak. 

Determining the resulting cosmogenic-activation back¬ 
ground contribution to a Oz//3/3-decay experiment re¬ 
quires estimating the production rates of the radioiso¬ 
topes in tellurium. Activation cross sections that span a 
wide range of neutron energies, from thermal up to sev¬ 
eral GeV, are therefore needed; however, experimentally- 
measured cross-section data is currently sparse. For neu¬ 
tron energies above 800 MeV, cross sections for neutron 
activation are expected to be approximately equal to 
those for proton activation, and can be estimated from 
existing experimental data for proton energies 800 MeV 
- 23 GeV [I7W5] . In these proton measurements, two 
long-lived radioisotopes were observed that have the po¬ 
tential to contribute background at the Oz//3/3-decay peak: 
ii° m Ag and 60 Co. Below 800 MeV, experimental data ex¬ 
ists for activation of natural tellurium by ~ 1-180 MeV 
neutrons 20 and activation of individual tellurium iso¬ 
topes by thermal to ~15 MeV neutrons ED; however, 
only a few reactions were measured, and no cross sections 
were reported for the production of 60 Co and 110m Ag. To 
deal with the lack of experimental data, the background 
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from cosmogenic activation has been estimated in the 
past (as in Ref. [7J) using a combination of the aforemen¬ 
tioned neutron and proton measurements and codes that 
either implement the semi-empirical formulae by Silber- 
berg and Tsao (S&T) [22H24] (e.g., YIELDX [22H24] . AC- 
TIVIA [25]), or are based on Monte Carlo (MC) methods 
(e.g., CEM03 [26], HMS-ALICE [2 7 \ , GEANT4 [28][29]). 

These estimates can be greatly improved with addi¬ 
tional neutron-activation cross-section measurements be¬ 
low 800 MeV, which can also be used to benchmark the 
S&T and MC codes. A sample of natural-Te 02 powder 
was irradiated at the Los Alamos Neutron Science Cen¬ 
ter (LANSCE) with a neutron beam containing neutrons 
with kinetic energies up to ^800 MeV, and having an en¬ 
ergy distribution that resembles the cosmic-ray neutron 
flux at sea-level. Following exposure, the 7 rays emit¬ 
ted from the sample were measured in a low-background 
environment with a high-purity-germanium (HPGe) de¬ 
tector to determine the radioisotopes present. Based on 
these results, flux-averaged cross sections were obtained 
for several dozen isotopes. 

The cross sections are used to investigate the impact 
cosmogenic activation will have on CUORE, a next- 
generation 0z^d/3-decay experiment that will use an ar¬ 
ray of 988 high-resolution, low-background natural-Te 02 
bolometers to search for the Oz^/3/5 decay of 130 Te. In 
addition, the measured cross sections are compared with 
cross sections calculated using the ACTIVIA code. De¬ 
tails of this measurement and subsequent analysis are 
discussed below. 



FIG. 1. Schematic of the target irradiated at LANSCE. The 
entire target is 6.2 cm long in the z direction. Each target 
component has cylindrical symmetry about the z-axis. This 
drawing is not to scale. Details on each component are given 
in Tabled] 


II. EXPERIMENTAL METHOD AND DATA 
ANALYSIS 

A. Target 

The target consisted of 272 g of natural-Te 02 powder 
held within a cylindrical plastic container wrapped on 
all sides with 0.05 cm of cadmium to remove thermal 
neutrons. The front and back cadmium-layers were also 
used to monitor the neutron flux on either side of the 
target. Circular aluminum and gold foils were placed 
throughout the target to monitor the neutron flux as well. 
The target geometry is illustrated in Figure |TJ and the 
details of each target component are listed in Table [I] 
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FIG. 2. LANSCE 30R neutron flux (red) [30] compared with 
the measured sea-level cosmic-ray-neutron flux (black) [31]. 


B. Neutron irradiation 

The target was irradiated with neutrons from the LAN¬ 
SCE 30R beam line for 43 hours during February 25-27, 
2012 . At LANSCE, neutrons are generated from spal¬ 
lation reactions induced by an 800 MeV pulsed proton 
beam incident on a tungsten target. The 30R beam 
line, which is 30° to the right of the proton beam, has 
a neutron-energy spectrum that closely resembles the 


cosmic-ray neutron spectrum at sea-level, but has an in¬ 
tensity 3 x 10 8 times larger, as shown in Figure^ A beam 
collimation width of 8.26 cm was used, which resulted in 
a beam-spot diameter of 8.41 cm at the target. 

The proton beam used to generate the neutrons con¬ 
sisted of 625-/is-long macropulses occuring at a rate of 
40 Hz. Each macropulse contained micropulses spaced 
1.8 fi s apart. The neutron time-of-flight was obtained 
by measuring the time between the arrival of the proton 
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TABLE I. Description of the target components illustrated in Figure [I] The material, dimensions, mass, and purpose of each 
component are given. The parameters A z and d are the thickness of the component along the z-axis and the diameter in the 
x-y plane, respectively. 


Component 

Material 

Az 

(cm) 

d 

(cm) 

Mass 

(g) 

Purpose 

Te0 2 

Te0 2 powder 

2.79 

6.43 

271.56 

Target 

All 

A1 

0.0813 

6.22 

6.68 

Neutron-flux monitor 

A12 

A1 

0.0813 

5.93 

6.06 

Neutron-flux monitor 

A13 

A1 

0.0813 

5.93 

6.06 

Neutron-flux monitor 

Aul 

Au 

0.00515 

2.54 

0.504 

Neutron-flux monitor 

Au 2 

Au 

0.00512 

2.54 

0.500 

Neutron-flux monitor 

Cdl 

Cd 

0.05 

6.7 

16.3 

Neutron-flux monitor 
Thermal-neutron absorber 

Cd 2 

Cd 

0.05 

7.3 

19.9 

Neutron-flux monitor 
Thermal-neutron absorber 

Cd 

Cd 

0.05 

— 

— 

Thermal-neutron absorber 

Plastic Container 

Polystyrene 

0.2 

— 

— 

Target holder 


macropulse at the tungsten target and the generation of 
a fission signal in a 238 U-fission ionization chamber [32] 
located 25.4 cm upstream of the TeC >2 target. The ioniza¬ 
tion chamber was only able to detect neutrons with en¬ 
ergies above the 238 U-fission threshold, which is approx¬ 
imately 1.25 MeV. The average neutron flux above 1.25 
MeV at the Te 02 target was determined to be 1.41 x 10 6 
neutrons/(cm 2 -s), with an estimated uncertainty of 10 % 
[33] based on uncertainties in the geometry and efficiency 
of the ionization chamber. 


C. Gamma-ray analysis of the irradiated target 

Approximately one week after the neutron irradiation, 
the Te 02 target was dismantled, and each component 
was analyzed using 7 -ray spectroscopy at the Lawrence 
Berkeley National Laboratory Low Background Facil¬ 
ity [34, |35]. The TeC >2 powder, cadmium foils, and 
aluminum foils were measured using an upright, 115%- 
relative-efficiency, n-type HPGe detector, and the gold 
foils were measured with a horizontal, 80%-relative- 
efficiency, p-type HPGe detector. Each detector was sur¬ 
rounded by a copper inner shield encased in a lead outer 
shield. The gold foils were highly activated and could be 
counted at a distance of 12 cm from the detector. The 
cadmium and aluminum foils had low levels of activity 
and were therefore measured directly on top of the de¬ 
tector to maximize the detection efficiency. For the Te 02 
powder, the 7 -ray measurements needed to be highly sen¬ 
sitive to long-lived radioisotopes, which had low levels of 
activity inside the powder. To maximize the detection 
efficiency, the Te 02 powder was mixed thoroughly and 
counted in a Marinelli beaker positioned over the top of 
the detector (Figure [ 3 ]). A plastic insert was placed in¬ 
side the beaker to decrease the thickness and increase 
the height of the powder, which in turn increased the 
solid angle of the detector seen by the powder and de¬ 
creased the self-attenuation of 7 rays from decays within 


the powder. The thickness and average height of the 
TeC>2 powder were 3.8 mm and ^5.6 cm, respectively. 
The TeC>2 was counted in this configuration periodically 
for six months to enable the observation of long-lived ac¬ 
tivation products after the short-lived ones decayed away. 
Figure [ 4 ] shows a 7 -ray spectrum for the TeC>2 powder 
collected four months after the irradiation. 


~5.6 cm 



* 

Detector 

endcap 


FIG. 3. Setup used during the 7 -ray measurement of the 
Te 02 powder. Each component has cylindrical symmetry 
about the dashed line. This drawing is not to scale. 


Each peak in the 7 -ray spectra was fit with a Gaussian 
summed with a quadratic background function to deter¬ 
mine the energy and net counts. For peaks with higher 
intensity, a smoothed step function was also added to the 
fitting function. The 7 -ray energies were used to identify 
the radioisotopes produced in the Te 02 powder. For 7 - 
ray lines that could come from the decay of more than 
one isotope, the contributors were identified from the de¬ 
cay half-life of the line. 

A list of the radioisotopes observed in the TeC >2 powder 
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FIG. 4. 3-day-long 7 -ray spectrum collected for the Te 02 
powder four months after the neutron irradiation, (a) Full 
spectrum, (b) A region of the spectrum where 110m Ag peaks 
were observed. Labeled peaks are associated with the decay 
of isotopes with Q values greater than the Q/ 3/3 of 130 Te, i.e., 
ii° m Ag (red and bolded), 124 Sb (blue), and 126 Sb (black). 

Other peaks in the region are from the decays of 125 Sb, 
129m Te? 105 Ag? and H4m In> 


is provided in Table \H\ Since 7 -ray measurements started 
one week after the neutron irradiation ended, only acti¬ 
vated isotopes with half-lives greater than day re¬ 
mained. Therefore, any observed isotope with a shorter 
half-life was a decay daughter of a longer-lived isotope. 
For example, the presence of 127 Te (9.35-hour half-life) 
and 129 Te (69.6-minute half-life) was due to the decays 
of the longer-lived metastable states 127m Te and 129 m Te, 
respectively. 


D. Photopeak efficiencies 

The 7 -ray measurements of the TeC>2 powder needed 
to be highly sensitive to long-lived radioisotopes, which 
had low levels of activity inside the powder. To maximize 
the detection efficiency, the powder was counted imme¬ 
diately next to the detector (Figure [ 3 ]). Determination 
of the photopeak efficiencies for the Te 02 powder from 
calibration measurements alone was impractical due to 
the complexity of the counting geometry and the effects 
of true-coincidence summing, which can be significant at 


such close range. Therefore, the efficiencies were obtained 
by running simulations with the Geometry and Tracking 
4 (GEANT4) code, version 4.9.4.p02, which were bench- 
marked against experimental measureme nts of various 
point and extended 7 -ray sources (Table III) that cov¬ 
ered a wide range of 7 -ray energies. 

For the benchmarking measurements, the 57 Co and 
54 Mn point sources were each counted at the center of 
the detector face and at four positions along the side of 
the detector that were spaced 2 cm apart and spanned 
the length of the HPGe crystal. The uranium source 
was counted on the side of the detector as well. Fol¬ 
lowing the natural-source method [36], the two extended 
sources, ESI and ES 2 , were constructed from powders 
that contained elements with naturally-occurring long- 
lived radioisotopes. ESI was designed to mimic the ge¬ 
ometry of the irradiated TeC>2 powder during the 7 -ray 
measurements, and ES 2 was designed to mimic both the 
geometry and density of the powder. Photopeak efficien¬ 


cies were obtained for all the 7 rays listed in Table III 


In addition, the total efficiency, which is needed to de¬ 
termine summing corrections, was obtained for the two 
57 Co 7 rays (122.06 keV and 136.47 keV) and the 54 Mn 
7 ray (834.85-keV). 

The benchmarking measurements were simulated using 
GEANT4. Each simulation included the HPGe detector, 
the 7 -ray source, and the lead and copper shielding. For 
each 7 ray of interest, the entire decay scheme of the par¬ 
ent nucleus was simulated. Angular correlations between 
coincident 7 rays were not taken into account; however 
at close distances to the detector, the effects on the pho¬ 
topeak efficiencies are largely averaged out and are thus 
small. 

Each simulated photopeak or total efficiency (e S7 ) was 
compared with the measured value (e m7 ), and the per¬ 
cent difference was determined: 


A = ernj e S7 x 100% _ 
G7 


( 1 ) 


Using the manufacturer’s detector specifications in the 
simulations resulted in Ae 7 values that ranged from ap¬ 
proximately -10% to -35%, with the agreement between 
simulation and measurement worsening at lower 7 -ray en¬ 
ergies. This kind of disagreement, especially overestima¬ 
tion by the simulation, has been seen in other studies that 
model the 7 -ray efficiencies of HPGe detectors using the 
geometry provided by the manufacturer (e.g., Refs. m 
SOI). Typically, the discrepancies have been attributed to 
physical characteristics of the detector (crystal location, 
Li-diffused-contact thickness, etc.) that are difficult for 
the manufacturer to precisely specify. When the source 
is counted close to the detector, small uncertainties in 
the detector’s parameters can have significant effects on 
the 7 -ray efficiencies. 


The adjustments listed in Table |TVj were applied to the 
detector geometry in GEANT4 to make the efficiencies 
from the simulations more closely match those from the 
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TABLE II. Radioisotopes observed in the irradiated TeCU powder. Unless otherwise indicated, all isotopes were produced by 
neutron interactions with tellurium. The measured and calculated flux-averaged cross sections (<73or and <ts&t, respectively) 
for neutron activation of tellurium are provided for isotopes with half-lives greater than 1 day. All ( 730 R were measured at the 
68% C.L. and are independent cross sections, except for those followed by “(cu),” which are cumulative. All rows corresponding 
to isotopes that can contribute background at the 130 Te 0^/3/3-decay peak are bolded, and for these isotopes, the decay modes 
(e and f3~ for electron capture and beta-minus decay, respectively) and Q values are given. 


Isotope 

Half-life 

<730R (68% C.L.) 

(mb) 

ffS&T 

(mb) 

Decay Q value 
(MeV) 

126 -pn 

12.93 d 




131 la| 

8.025 d 




lls Te 

6.00 d 

5.7 zb 1.2 

9.80 


119mrj-^ 

4.7 d 

6.3 ± 0.8 

13.0 

2.554 (e) 

121 rpJb I 

19.17 d 




121 

164.2 d 

16 ±2 

25.2 


123mrp^ 

119.2 d 

36 ±4 

8.6 


125mrp^ 

57.4 d 

83 ± 10 

17.0 


127 Te 

9.35 h 




127mrp^ 

106.1 d 

46 ±9 

25.3 


129 Te 

69.6 m 




129mrp^ 

33.6 d 

53 ± 17 (cu) 

22.4 (cu) 


131 Te 

25 m 




131mrpjb] 

33.25 h 




131m-yr Jal 

11.84 d 




118 Sb 

3.6 m 



3.657 (e) 

119 sfl 

38.19 h 




120 4 

5.76 d 

6.3 ± 0.8 

10.2 

2.681 + E ex (e) 

122 Sb 

2.7238 d 

14 ± 2 (cu) 

15.4 (cu) 


124 Sb 

60.2 d 

16 ± 2 (cu) 

19.1 (cu) 

2.904 (/3~) 

125 Sb 

2.759 y 

18 ± 2 (cu) 

18.8 (cu) 


126 Sb 

12.35 d 

6.7 ± 0.9 (cu) 

26.4 (cu) 

3.673 (f3~) 

127 Sb 

3.85 d 

13 ± 2 (cu) 

9.8 (cu) 


113 Sn 

115.1 d 

2.6 ± 0.3 (cu) 

3.0 (cu) 


117m Sn 

14 d 

4.3 ±0.6 

0.63 


111 In 

2.805 d 

2.3 ± 0.3 (cu) 

2.1 (cu) 


114m In 

49.51 d 

1.9 ±0.2 

0.31 


105 Ag 

41.29 d 

0.56 ± 0.07 (cu) 

0.45 (cu) 


106m^g 

8.28 d 

0.44 ± 0.09 

0.39 

3.055 (e) 

110 Ag 

24.56 s 



2.893 (/3~) 

110m Ag 

249.83 d 

0.28 ± 0.04 

0.054 

3.010 (f3~) 

ln Ag 

7.45 d 

0.42 ± 0.09 (cu) 

0.030 (cu) 


101 Rh 

3.3 y 

0.06 ± 0.01 (cu) 

0.24 (cu) 


101m Rh 

4.34 d 

0.30 ± 0.05 (cu) 

0.24 (cu) 


102m Rh 

3.742 y 

0.15 ±0.02 

0.12 


60 Ccpl 

5.27 y 

< 0.0016 (cu) 

0.0013 (cu) 

2.823 (/3~) 

7B fr 

53.24 d 

1.4 ±0.2 

2.5 



a This isotope was produced by interactions with spallation protons created in the target during the neutron irradiation. Therefore, no 
cross sections are provided. 

b This isotope had a high probability of being produced by interactions with < 1.25-MeV neutrons. Therefore, no cross sections are 
given. 

c A flux-averaged cross section could not be obtained for 119 Sb because the strongest 7 -ray line at 24 keV overlapped with x-rays 
emitted by other activated isotopes. 

d 60 Co was not conclusively observed in the 7 -ray spectra due to 102m Rh and 110m Ag peaks being present where the 60 Co peaks were 
expected. Therefore the cross section quoted for 60 Co is an upper limit. 

e 7 Be was produced almost exclusively by neutron interactions with oxygen. The cross sections given correspond to these interactions. 


benchmarking measurements. The larger disagreement 
at low energies between the simulated and measured ef¬ 
ficiencies pointed to additional, unspecified attenuating 
material that was present in the actual detector. To ad¬ 


dress this, the thickness of the aluminum mounting cup 
that immediately surrounds the HPGe crystal was in¬ 
creased by 2.25 mm to achieve closer agreement between 
the simulations and measurements. 
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TABLE III. Description of 7 -ray sources used to benchmark GEANT4. 


Source 

Composition 

Dimensions 

7 -ray 

Branching Ratio 




(keV) 

(%) 

co-57 

co-57 

Point source 

122.06 

85.60 ±0.17 




136.47 

10.68 ±0.08 

Mn-54 

Mn-54 

Point source 

834.85 

99.9760 ±0.0010 

Uraniunj^] 

Natural uranium ore (0.1176 g) 

Diameter = 4.76 cm 

185.72 ( 235 U) 

57.2 ±0.8 

mixed with epoxy 

Thickness — 3.175 mm 

46.54 ( 210 Pb) 

4.25 ±0.04 




186.21 ( 226 Ra) 

3.64 ±0.04 




242.00 ( 214 p b) 

7.251 ±0.016 




295.22 ( 214 p b) 

18.42 ± 0.04 




1764.49 ( 214 Bi) 

15.30 ±0.03 

ES|] 



2204.06 ( 214Bi ) 

4.924 ±0.018 

La 203 powder (89 g), 

Inner radius = 5.06 cm 

201.83 ( 176 Lu) 

78.0 ±2.5 

LU 2 O 3 powder (2 g), 

Outer radius = 5.443 cm 

306.78 ( 176 Lu) 

93.6 ± 1.7 


KC1 powder (4 g) 

Average height = 5.75 cm 

788.74 ( 13 S La) 

34.4 ±0.5 




1435.80 ( 138 La) 

65.6 ±0.5 




1460.82 ( 40 K) 

10.66 ±0.18 




269.46 ( 223 Ra) 

13.9 ±0.3 




271.23 ( 219 Rn) 

10.8 ± 0.6 




832.01 ( 211 p b) 

3.52 ±0.06 

ESP 



351.07 ( 211 Bi) 

13.02 ±0.12 

(Unirradiated) Te 02 powder (228 g), 

Inner radius = 5.06 cm 

Note: All 7 rays used 



La 203 powder (23 g), 

Outer radius = 5.443 cm 

to analyze ESI were also 



LU 2 O 3 powder (6 g), 

K 2 SO 4 powder (14 g) 

Average height = 6.5 cm 

used to analyze ES2. 



a All isotopes in the source were assumed to be in secular equilibrium. 

b Due to a small 227 Ac contamination in the La 2 C> 3 , ESI and ES 2 also contained 227 Ac and its daughter isotopes, which were assumed 
to be in secular equilibrium with each other. 7 -rays from the 227 Ac chain are also listed in the table. 


Figure [5] shows the values of Ae 7 obtained after the 
adjustments to the detector geometry were made in the 
GEANT4 simulations. The uncertainties in Ae 7 take into 
account the statistical uncertainties in the measurements 
and the simulations, as well as the uncertainties in the 
source activities and branching ratios of the 7 rays. The 
total uncertainty in the simulated efficiencies was esti¬ 
mated to be 5%, which is slighly larger than the standard 
deviation of Ae 7 . 

The photopeak efficiencies of the 7 rays used to identify 
the isotopes in Table [II] were obtained for the irradiated 
Te 02 powder by performing GEANT4 simulations using 
the adjusted detector values in Table [IVJ Simulations in¬ 
dicate that summing could have as much as a 40% effect 
for certain photopeak efficiencies. Figure [5] gives confi¬ 
dence that the GEANT4 simulations could model sum¬ 
ming correctly and provide photopeak efficiencies for the 
irradiated TeC >2 powder with around 5% uncertainty. 


E. Neutron transmission during neutron 
irradiation 


Following the neutron irradiation, the aluminum and 
cadmium foils located in front of and behind the Te 02 
powder were measured with an HPGe detector, as was de¬ 


scribed in Section IIC The total neutron transmission, 



FIG. 5. Percent differences between the measured and sim¬ 
ulated 7 -ray efficiencies as a function of 7 -ray energy. The 
simulated efficiencies were obtained using the adjusted val¬ 
ues in Table |TV[ Points corresponding to total efficiencies are 
indicated with “(Total)” in the legend. All other points cor¬ 
respond to peak efficiencies. 


T toU through the Te 02 powder was estimated by compar¬ 
ing the activities of the activation products in the front 
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TABLE IV. Detector parameters adjusted in the GEANT4 simulations. The nominal values provided by the manufacturer are 
given, along with the values that allowed for satisfactory (~5%) agreement between the efficiencies from the simulations and 
the benchmarking measurements. 


Parameter 

Nominal Value 

Adjusted Value 


(mm) 

(mm) 

Length of HPGe crystal 

85.5 

80.f\ 

Distance between HPGe crystal and detector window 

0 


Thickness of aluminum mounting-cup 

0.5 

2.75 

Thickness of internal dead layer (lithium contact) 

1 

2 


a The effects of shortening the crystal in the simulation could also be reproduced by using the nominal length of the crystal and adding 
a 1.85-mm-thick, 3.5-cm-long copper ring around the aluminum mounting-cup, 8 cm below the top of the detector endcap. Since the 
presence of such a ring was not specified by the manufacturer and the corresponding simulations provided equivalent results to the 
shortened crystal geometry, the simulated efficiencies for the 80.5-mm-long crystal were used in the cross-section analysis. 


TABLE V. Neutron-transmission results. Reactions used to 
determine the neutron transmission through the TeC >2 powder 
are listed, along with the corresponding values of T to t- 


Reaction 

Ttot 

27 Al(n,X) 22 Na 

0.98 ±0.03 

Cd(n,X) 105 Ag 

0.84 ±0.01 

Cd(n,X) 110m Ag 

0.86 ±0.01 


foils (All and Cdl) with the activities in the back foils 
(A13 and Cd 2 ). The values of T tot obtained by analyz¬ 
ing the activation products are given in Table \V\ and the 
uncertainties quoted are statistical. To account for the 
variation in the results, T tot was taken to be 0.90 ± 0 . 10 , 
which spans the range of values given in Table [V] along 
with their uncertainties. The average neutron transmis¬ 
sion, T, through the Te 02 powder was then estimated to 
be 


T ~ 1 + — = 0.95 ± 0.05. 
2 


( 2 ) 


B 7 is the branching ratio of the 7 ray, and e 7 is the pho¬ 
topeak efficiency of detecting the 7 ray. The production 
rates for 125m Te and 127m Te were described by more com¬ 
plex growth-and-decay relations and were obtained using 
the appropriate modifications to Equation [3] 


G. Flux-averaged cross sections 

The flux-averaged cross section, 03011 , for neutron ac¬ 
tivating an isotope in the irradiated Te 02 powder is de¬ 
termined from 


f%™<T(E)'p 30K (E)dE 

/|T>30R(£)d£ 


(4) 


where cr(E) is the cross section for producing the iso¬ 
tope with neutrons of kinetic energy E , (P 3 or(E) is the 
differential neutron flux hitting the front of the target 
in units of [neutrons/(cm 2 -s-MeV)], and E m [ n and E max 
are respectively the lowest and highest neutron energies 
hitting the TeC >2 powder. 

The isotope-production rate can also be expressed as 


F. Isotope-production rates 

The production rate for each isotope can be deter¬ 
mined using data from the 7 -ray spectra collected for the 
irradiated TeC >2 powder. In most cases, the isotope pro¬ 
duced in the powder is not fed by other isotopes during 
or after the neutron irradiation. Under this condition, 
the production rate, i? 30 R, can be obtained using: 

#30R = g J ? (3) 

where A is the decay constant of the isotope, C 7 is the 
number of counts in the 7 -ray peak of interest corrected 
for the growth and decay of the isotope during the irra¬ 
diation and the decay of the isotope after the irradiation, 


pE max 

R 3 or ~ N&30RT / (f30R(E)dE , ( 5 ) 

J -Emin 

where N is the number of tellurium nuclei in the powder 
(except for the production of 7 Be, where N is the number 
of oxygen nuclei in the powder). 

The total neutron flux below 1.25 MeV during the ir¬ 
radiation was determined to be nonnegligible from the 
amount of 198 Au created by (n/y) reactions in the gold 
foils. Therefore, cross sections could be obtained only 
for isotopes produced soley (or primarily) by interactions 
with neutrons of energy > 1.25 MeV, and in these cases, 
Emin and F7 max from Equation [5] could be set to 1.25 MeV 
and 800 MeV, respectively. 

The flux-averaged cross sections, shown in Table [ilj 
can then be determined from Equations [3] and [5] 















III. COMPARING MEASURED AND 
CALCULATED CROSS SECTIONS 

Isotope-production cross sections for tellurium were 
also obtained by using the ACTIVIA code to perform 
calculations based on the S&T semi-empirical formu¬ 
lae. These formulae were originally developed to describe 
proton-nucleus interactions, but they are assumed to be 
applicable to neutron-nucleus interactions as well. The 
calculated cross sections are reported in Table [ilj Al¬ 
though the formulae are only valid for proton and neu¬ 
tron energies > 100 MeV and they do not distinguish 
between ground and metastable states in product nuclei, 
the calculated and measured cross sections agree reason¬ 
ably well, within a factor of 3 on average. One should 
note that the cross section calculated for 110m Ag was un¬ 
derestimated by approximately a factor of 5. 

IV. COSMOGENIC-ACTIVATION 
BACKGROUND IN THE CUORE EXPERIMENT 

The CUORE experiment will use an array of 988 high- 
resolution, low-background Te 02 bolometers to search 
for the Ovpp decay of 130 Te. Each bolometer is com¬ 
prised of a 5x5x5 cm 3 natural-Te 02 crystal that serves 
as both a source and a detector of the decay. CUORE is 
aiming for a background rate of 10 -2 counts/(keV-kg- y) 
at the 130 Te value of 2528 keV, which would allow 
the experiment to reach a half-life sensitivity of 9.5 x 10 25 
years (90% C.L.), assuming a live time of 5 years and a 
full-width-at-half-maximum energy resolution of 5 keV 
011 - 

Using the results of the neutron-activation measure¬ 
ment discussed in this work and the proton-activation 
measurements of Ref. [19], one can determine the back¬ 
ground contribution to CUORE from the cosmogenic ac¬ 
tivation of the Te 02 crystals that occurs during sea trans¬ 
portation from the crystal-production site in Shanghai, 
China to LNGS in Italy. The results of both this work 
and Ref. m indicate that 110m Ag and 60 Co are the only 
two long-lived radioisotopes that will contribute mean¬ 
ingfully to the background at the Ou(3 /3-decay peak due 
to their Q values being greater than Q/ 3 / 3 . 110 Ag will also 
contribute a small amount to the background because 
110m Ag decays to it 1.33% of the time. 

The production rates, R, of 110m Ag and 60 Co were each 
estimated to be 


parameter F B syd from Ref. [31] taken to be 0.73 =b 0.22 
[33] for the route used to ship the TeCU crystals. 80% 
of the 110m Ag and as much as 37% of the 60 Co were 
produced by 1.25-800 MeV neutrons. 

The fraction of 110m Ag, 110 Ag, and 60 Co decays that 
deposit energy in a 60 keV-wide region-of-interest (ROI) 
surrounding the 0zqd/3-decay peak was estimated using 
GEANT4 simulations of a single 5x5x5 cm 3 TeCU crys¬ 
tal. The values obtained were 0.5%, 0.4%, and 1% for 
ii° m Ag, 110 Ag, and 60 Co decays, respectively. In the 
full CUORE array, the presence of nearby crystals would 
often lead to energy being deposited in more than one 
crystal. As most 0zqd/%decays would deposit all of their 
energy in a single crystal, the background can be reduced 
by rejecting events in which energy was deposited in more 
than one crystal. Simulations of a 3 x 3 x 3 Te02-crystal 
array indicate that rejecting multi-crystal events can sup¬ 
press the 110m Ag contribution to the ROI by a factor of 
~2, while the contributions from 110 Ag and 60 Co will be 
minimally affected. 

To estimate the background rate at Qpp from cosmo¬ 
genic activation of Te 02 , the following assumptions were 
made: (1) each crystal spends 3 months at sea level, (2) 
no 110m Ag, 110 Ag, and 60 Co were present at the begin¬ 
ning of shipment due to their removal during the crystal- 
growth process, and (3) crystals were delivered to LNGS 
and stored underground at a constant rate from early 
2009 to late 2013 [33]. The resulting contamination levels 
for 110m Ag+ 110 Ag and 60 Co when CUORE begins oper¬ 
ation in late 2015 will be ^2 x 10 -8 Bq/kg and ^10 -9 
Bq/kg, respectively, which correspond to background 
rates of ^6 x 10 -5 counts/(keV-kg- y) and x 10 -6 
counts/(keV-kg• y), respectively. After 5 years of run¬ 
ning, the contamination levels will decrease to ^2 x 10 -1 ° 
Bq/kg for 110m Ag-f- 110 Ag and ^6 x 10 -1 ° Bq/kg for 
60 Co, which correspond to background rates of ^4 x 10 -7 
counts/(keV-kg• y) and ^4x 10 -6 counts/(keV-kg- y), re¬ 
spectively. The contamination levels given here are lower 
than those predicted in Ref. [7] due to Lozza et al. as¬ 
suming a longer exposure time of 1 year and a shorter 
overall cooling time underground of 2 years. Reject¬ 
ing multi-site events should decrease the 110m Ag+ 110 Ag 
background rates by a factor of ^2. Although the back¬ 
ground rates in the ROI are at least two orders of mag¬ 
nitude lower than the current CUORE goal background 
of 10 -2 counts/(keV-kg- y), for future experiments striv¬ 
ing for essentially zero background, cosmogenic activa¬ 
tion may have to be addressed more stringently. 


CR,i, (6) 

i 


V. CONCLUSIONS 


where cq is the isotope-production cross section assigned 
to energy bin i, and (f)cR,i is the differential cosmic-ray 
neutron flux at sea-level integrated over energy bin i. The 
energy bins, integrated fluxes, and cq values are given in 
Table VI The cosmic-ray neutron flux determined by 
Gordon et al. ED was used in this analysis, with the 


Flux-averaged cross-sections for cosmogenic-neutron 
activation of radioisotopes in natural tellurium were mea¬ 
sured by irradiating Te 02 powder with a neutron beam 
containing neutrons of kinetic energies up to ^800 MeV, 
and having an energy spectrum similar to that of cosmic- 
ray neutrons at sea-level. The cross sections obtained for 
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TABLE VI. Energy bins used in the estimation of the 110m Ag and 60 Co production rates. The differential cosmic-ray neutron 
flux at sea-level integrated over each bin is provided. The isotope-production cross sections assigned to each bin are also listed. 
For bin 1, the cross sections obtained in this work are used. For bins 2, 3, and 4, the cross sections used were those measured 
in proton-activation experiments with 800 MeV, 1.4 GeV, and 23 GeV protons respectively. The individual contributions to R 
in units of [s 1 ] and [%] are given in the last two columns. 


Bin 

Bin Range 

Integrated Neutron Flux 
(s 1 -cm 2 ) 

Cross section 
(mb) 


Contribution to R 

(s- 1 4 ) 






110m Ag 

60 Co 


110m Ag 


60 Co 

1 

1.25 MeV 

- 800 MeV 

(3.7 ± 1.3) 

x 10 -3 

0.28 ±0.04 

< 0.0016 

(2.9 ± 1.1) x 

10 ' 

- 6 < (1.7 ±0.6) x 10“ 8 * 









(80%) 


(< 37%) 

2 

800 MeV 

- 1.4 GeV 

(5.3 ± 1.9) 

x 10 -5 

3.95 =b 0.4CF1 02] 

0.09 ±0.04 

m 

(5.9 ±2.2) x 

10 " 

_7 (1.4 ± 0.8) x 10 -8 






1—1 



(16%) 


(> 30%) 

3 

1.4 GeV 

- 23 GeV 

(2.6 ±1.0) 

x 10“ 5 

1.9 ±0.3 ns] 

0.20 ±0.04 

nn 

(1.4 ±0.6) x 

10 ' 

' 7 (1.5 ± 0.6) x 10 -8 









(3.9%) 


(> 33%) 

4 

23 GeV - 

150 GeV 

(1.6 ±0.6) 

x 10“ 7 

0.88 ± 0.59 pi] 

0.75 ±0.08 

m 

(4.0 ±3.1) x 

10 “ 

10 (3.4 ± 1.3) x 10“ 10 









(0.01%) 


(> 0.8%) 


a The value of this cross section was reported incorrectly in Ref. ns, but correctly in Ref. [42]. 


nomAg an d 60 Co, the two isotopes which have both half- 
lives of order a year or longer and Q values larger than the 
Qpp of 130 Te, were combined with results from tellurium 
activation measurements with 800 MeV - 23 GeV protons 
to estimate the background in the CUORE experiment 
from cosmogenic activation of the Te02 crystals. The 
anticipated 110m Ag+ 110 Ag and 60 Co background rates in 
[counts/(keV-kg- y)] at the O^/3/3-decay peak were deter¬ 
mined to be ^6 x 10 -5 and ~7 x 10 -6 , respectively, at the 
beginning of counting and ^4 x 10 -7 and ^4 x 10 -6 , re¬ 
spectively, after 5 years of counting. The 110m Ag+ 110 Ag 
rates should decrease by a factor of ^2 if multi-crystal 
events are efficiently rejected. These rates are at least 
two orders of magnitude lower than the goal background 
for the CUORE experiment. 
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